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C^ . Abstract 



We investigate the crumpling transition for a dynamically triangulated random surface embedded 

in two dimensions using an effective model in which the disordering effect of the X variables on the 

pH \ correlations of the normals is replaced by a long-range "antiferromagnetic" term. We compare the 

results from a Monte Carlo simulation with those obtained for the standard action which retains the 

X's and discuss the nature of the phase transition. 



It is now clear that the pathologies of the continuum Polyakov action for the bosonic string |l|] in 
physical dimensions (d > 1) 

Z^ f DgDX exp (-^f d^x^g'^'daX^d^X^^ (1) 

are also manifest in the discretized Euclidean form of the action used in simulations 

^ = E /n^^^-pR E i^t-xf;A . (2) 

T •' i \ <ij> J 

In this simple gaussian action the continuum integration over the metric is replaced by a summation over 
triangulations of the discretized worldsheet, the variables Xi which specify the embedding of the surface 
typically live at the nodes of the triangulation and the sum is over the edges < ij >. Both microcanonical 
simulations, in which the number of nodes was fixed, and grand canonical simulations, in which it varied, 
were carried out and it was discovered that the worldsheets generated by the action were very crumpled 
and failed to give a non-trivial continuum limit H. 

A possible remedy, suggested by Polyakov in H, was the inclusion of a stiffness term in the action 
which could be written using either of the two equivalent forms of the extrinsic curvature squared 

K^ = f d'^xVh{AhX''f 

= f d^xVhh''''VanfVbnl,. (3) 

Subsequent simulations revealed that discretizations of the first version were also afflicted by pathologies 
Q but that a discretization of the second form 

Se^ Y. (l-"A,.nA,), (4) 

where the A's are adjacent triangles and the n's normals, offered the possibility of defining a continuum 
limit at a second order phase transition where the string tension scaled [pj . Some recent simulations have 
cast doubt on the second order nature of the transition, but it appears that the requisite scaling of the 
string tension may still be taking place M . Other alternative actions are also being investigated M . 
The extrinsic curvature squared term in eq.H may also be written as 

Se^ Y. (1-C0S%), (5) 

where the 9ij is the angle between the adjacent normals. In the degenerate case of a surface embedded 
in just two dimensions it is still possible to write this term using Ising-like variables a = ±1 

Se= Y (l-^^'^j) (6) 

<v> 

as the cosine will always take the value ±1. (Note that these Ising variables live on the dual of the 
triangulation.) This prompted Renken and Kogut to suggest that if one ignored the effect of the X 
variables in the gaussian term one might expect an Ising like crumpling transition in this case as 5*6 
in eq.H is essentially the 2d Ising model action. This transition would thus be second order and might 
serve as a model for the transition in more realistic dimensions. However, their numerical simulations of a 
gaussian plus extrinsic curvature squared action on a fixed triangulation using finite size scaling techniques 
revealed a first order transition S. It does appear that this model on a dynamical triangulation has a 
second order crumpling transition which is more in line with expectation, though it is only possible to 
work with much smaller meshes in this case [Q . 

As we have ignored the gaussian term it is natural to enquire what effect its inclusion might have on 
the analysis above. It was pointed out in |1^ that if one used a proper time representation for the X 
propagator of a gaussian surface action 

<X^{i^)X'^{^,)>^6^'^ r i^cxp(-a(e.-C,)^), (7) 



where L is an infrared cutoff and A an ultraviolet cutoff, then the normal-normal correlations were of the 
form 

<^te)-'^(0)>- -(^-ry^- (8) 

The constant A depends on the dimension in which the surface is embedded [J. From eq.|8| it is clear 
that the gaussian term induces an antiferromagnetic interaction between the normals. This suggests a 
possible interpretation of the crumpling transition as arising from the competition between the gaussian 
"antiferromagnetic" term and the extrinsic curvature "ferromagnetic" term, which tends to order the 
normals. 

In this paper we test this interpretation for a surface embedded in two dimensions by comparing the 
behavior of the "effective" action where the gaussian term has been replaced by a term mimicing its effect 
on the normal spins, 

S = AY^^a,a,+\Y,{l-a.,a,), (9) 



r 



with the standard action we simulated in [pi 

where A is the ferromagnetic coupling proportional to the inverse temperature. For convenience we set 
the constant A = 1 in the simulations. Note that the sum in the first term of eq.0 is over all the spins 
so we have a long-range interaction. Moreover r^ is the intrinsic distance between the spins <7i and aj , 
i.e. the shortest path (number of links) between points i and j in the triangulation. In fact since there 
are no X variables in eq.|9|, the triangulation is no longer embedded in any dimension and there are no 
extrinsic distances. This means that we can no longer measure the radius of gyration to see how the size 
of the surface is changing, though we can measure the contribution of the extrinsic curvature squared 
term eqJq which we would expect to be small in a smooth phase and large in a crumpled phase. 

The action in eqJfl is rather similar to that suggested in pT| to explain the phenomenon of spin canting 
and reentrance g 

S^-JoJ2s^■ Sj +2eJ2 -^S. ■ Sj (11) 

where the S's are classical two-component spins and the sum in the second term is restricted to rij < R, 
where R is some cutoff. The authors in |ll] were interested in d = 3 and a < 2 and argued that the 
competition between the short range Jq term and the long-range e term accounted for canting. Actions 
of the form 

S = -Joy^^—-a,aj (12) 



ij 



r, 



u 



have also been considered [|2[, and display phase transitions for d < a < 2d. In our case the l/r'* 
interaction induced between all the normals is apparently on the borderline a = 2d if we assume d — 2, 
but direct measurement of d on the highly irregular meshes generated by 2d quantum gravity that we use 
gives d > 2. However our interaction is antiferromagnetic, so it appears that the analogy with the canting 
action is closer as we have two competing terms, one "ferromagnetic" and one "antiferromagnetic" , and 
observe a transition as we vary their relative strengths. 

The most difficult part in performing a Monte Carlo simulation of the action in eq.0 is calculating 
the intrinsic distance r^ (represented as a matrix) between all points i,j on the random triangulation. 
There is a straight-forward algorithm which does this but it is 0{N'^), where N is the number of points 
in the triangulation. This makes it prohibitively expensive to simulate dynamical random triangulations, 
since after every change to the triangulation one would have to recalculate r^ (even though the changes 
to the triangulation would consist only of the standard "flip" move which is local this would still produce 
non-local changes in the intrinsic distance matrix r^). Therefore all our simulations use fixed random 

^The considerations of | |lO[ | were for rigid or crystalline surfaces but they are essentially unchanged for dynamical surfaces. 

^ Reentrance is the behavior displayed by metallic spin glasses such as AuFe which have a paramagnetic to ferromagnetic 
transition at some T^ and a ferromagnetic to spin-glass transition at some Tf < T^. Below Tj spins tend to rotate away 
from an external field, which is termed canting. 



triangulations so that we have only to calculate rij once at the beginning. As our model is essentially an 
Ising model coupled to two-dimensional quantum gravity, albeit with an extra long-range term, we have 
chosen to use random triangulations coming from pure two-dimensional quantum gravity simulations as 
was done in |T^. 

We have performed simulations on three sizes of random triangulation of spherical topology with 
total number of points N — 100, 500 and 1000. For each N we ran at roughly 20 values of A between 1 
and 3. After thermalizing we do 100000 Monte Carlo updates using the standard Metropolis algorithm, 
measuring the energy, magnetization and spin-spin correlation function after every update which allows 
us to calculate autocorrelation times and correlation length. We find, as expected for the local Metropolis 
algorithm, that the autocorrelation time r scales as the square of the correlation length ^; in fact for the 
energy, by fitting r ~ ^^, we extract a dynamical critical exponent z = 2.1(2). From the fluctuations in 
the energy and magnetization we obtain the specific heat and susceptibility in the usual manner. 

In Fig. 1 we show the total energy E from the largest system simulated N — 1000 (graphs from 
A^ = 100 and 500 look almost identical) and its two parts - Eir is the antiferromagnetic long-range part 
(first term in eq.0 with A=l) and i?„„ is the ferromagnetic nearest-neighbor part (second term in eqM 
including A). In Fig. 2 we show the total magnetization Af , from all three simulations. At small A, which 
corresponds to high temperature, we see from Fig. 1 that most of the energy is in the ferromagnetic piece 
Emi since the spins are disordered, which we would expect for a "crumpled" phase - M is small in Fig. 
2. Conversely, at large A (low temperature) the spins line up so M ^ 1 and most of the energy resides 
in the antiferromagnetic Eir- In this smooth phase the contribution from the extrinsic curvature squared 
Enn is small. From these we deduce there is some sort of crumpling transition at around A = 2.3 — 2.4. 
Also from Fig. 2 we see that as the system size increases the "jump" in M becomes more abrupt, perhaps 
signaling the appearance of a first order phase transition. However there is no corresponding jump in the 
total energy so this seems unlikely. 

Stronger evidence of the phase transition not being of first order comes from the specific heat, shown 
in Fig. 3. For a first order phase transition, standard finite-size scaling theory predicts that for a 
two-dimensional system of size N = L x L, the specific heat peak 



whereas for a second order transition 



Cma. = AL^ + B, (13) 



Crr,a.=A'L"/''+B'. (14) 



We do not know a priori that our system is two-dimensional so we must write L — Nd , where d is 
the fractal dimension of the random triangulation. In fact, numerical simulations for pure 2d quantum 
gravity [Q and for 2d quantum gravity coupled to Potts models |1^ yield d ~ 2.7 or 2.8; and an analytical 
calculation predicts that d lies between 2 and 3 for pure 2d quantum gravity |1(| . Thus we fit the specific 
heat peak to 

Cma. = A'N'^/'^" + B' (15) 

for our three values of N obtaining a/vd = 0.12(2). Thus a first order phase transition is apparently 
ruled out. One could argue that our systems are too small to see the true scaling of Cmax, however 
Renken and Kogut managed to see its first order scaling in their model on systems of size 576 and 1024 
[p|. It should be noted, however, that they carried out their simulations on a regular, fixed lattice rather 
than the highly disordered lattice we are using, so the different behavior we observe might be due to the 
lattice rather than the absence of the gaussian term. 

If we assume that the phase transition is second order and use the scaling relation 

vd^2-a (16) 

we obtain a — 0.20(3) and vd ~ 1.8(3). Now, from our measurements of the correlation length, we can 
obtain v by fitting 

C'-IAc-AI-- (17) 

for A < Ac, with Ac = 2.35. This value for Ac is obtained from the peak in the susceptibility, shown in 
Fig. 4. This leads to i^ = 0.7(1). Therefore we estimate d — 2.6(8). Unfortunately the error here is 
rather large but nevertheless the numbers are consistent. As a final check we could assume a = (as is 



the case for the phase transition in the usual 2d Ising model) , then vd — 2 so we eliminate one source of 
error and obtain d = 2.9(4). 

At this point we could tentatively conclude that our data is consistent with the phase transition being 
second order. However, it is interesting to note similarities with data from numerical simulations of 
(gaussian) spin glasses which typically have a rather broad peak in the specific heat and a susceptibility 
peak at lower temperature |l^ (which corresponds to higher A). If we look at Fig. 3 and Fig. 4 we see just 
this behavior. The very weak growth of the peak in the specific heat that we have measured could be 
a finite size effect masking a negative a (in the Sherrington/Kirkpatrick spin glass model a = — 1 ]l7t) 
n. Moreover, at the phase transition in our model we have competing interactions which typically leads 
to spin-glass type behavior. As these interactions are still present in the original model with gaussian 
plus extrinsic curvature-squared terms this suggests that the low A, crumpled phase in this model could 
conceivably be spin-glass like rather than paramagnetic. A simulation by Heermann [nsl has, in fact, 
found that the radius of gyration X2 is non self-averaging in the crumpled phase, a behavior found in 
some observables in spin glasses. 

To summarize the numerical results and speculations: we have simulated an effective model for 
crumpling in two dimensions in which the X variables are replaced by an antiferromagnetic term that 
mimics their effect on the normals. We have found what appears to be a weak second (or possibly higher) 
order transition, rather than the first order transition of [||, which retained the Xs. We have pointed 
out that this might be due to the differing underlying fixed meshes in the simulations rather than the 
absence of the X variables. We have also speculated that the competing nature of the interactions at 
the transition may give a spin-glass like phase at low A rather than a paramagnetic phase and presented 
some numerical support for this view. 

It would be interesting to repeat the simulation on a fixed regular mesh to compare directly with 
|g[ and to extend the effective model work to the more realistic case of a surface embedded in three 
dimensions, where one would employ continuous spins. The nature of the low A ("crumpled") phase also 
merits further investigation for both the effective model and the original action. 

This work was supported in part by NATO collaborative research grant CRG910091. CFB is supported 
by DOE under contract DE-AC02-86ER40253 and by AFOSR Grant AFOSR-89-0422. 



^Recent very high statistics simulations of a gaussian plus extrinsic curvature-squared action in three dimensions [h| have 
also found very little, if any, growth in the peak. 
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Figure Captions 

Fig. 1. Energy for the N=1000 simulation. 

Fig. 2. Magnetization for all three (N=100,500 and 1000) simulations. 

Fig. 3. Specific heat for all three simulations. 

Fig. 4. Susceptibility for all three simulations. 



